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Trigonal warping and anisotropic band splitting in monolayer graphene due to Rashba 

spin-orbit coupling 
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We study the electronic band structure of monolayer graphene when Rashba spin-orbit coupling 
is present. We show that if the Rashba spin-orbit coupling is stronger than the intrinsic spin-orbit 
coupling, the low energy bands undergo trigonal-warping deformation and that for energies smaller 
than the Lifshitz energy, the Fermi circle breaks up into separate parts. The effect is very similar 
to what happens in bilayer graphene at low energies. We discuss the possible experimental implica- 
tions, such as threefold increase of the minimal conductivity for low electron densities, anisotropic, 
wavenumber dependent spin-splitting of the bands and the spin polarization structure. 
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PACS numbers: 73.22.Pr,73.23.Ad 

Follovifing the pioneering work of Kane and Mele^ 
where a novel type of quantum spin-Hall effect was pre- 
dicted, there has been a significant interest in studying 
the spin-orbit (SO) coupling in graphene. Generally, the 
SO interaction in graphene can have intrinsic or extrin- 
sic origins. Density functional and other calculations 
indicate^""— that the intrinsic spin orbit (ISO) coupling 
is in the range of 1 — 50 /ieV and therefore it could only 
be probed at very low temperatures. The SO interaction 
can also be induced externally, by impurities^, by ex- 
ternal electric field perpendicular to the graphene sheet 
or by electrostatic interaction with the substrate^. The 
latter two mechanism can introduce Rashba type spin- 
orbit (RSO) interaction, which will be in the focus of our 
interest in this paper. Many experimental implications 
of both intrinsic and extrinsic SO coupling in graphene 
have also been addressed in the literatureii^"— . The 
recent experiment of Varykhalov et al^ have shown that 
it was possible to fabricate a quasi freestanding graphene 
on Ni(lll) surface with a single layer of intercalated gold 
atoms between the graphene layer and the Ni surface. A 
spin-orbit interaction induced band splitting of « 13 meV 
was measured, which is two-three orders of magnitude 
larger than the predictions for ISO^"— . The spin-orbit 
coupling was identified as RSO and the large splitting 
was attributed to the high nuclear charge of the interca- 
lating gold atoms. Very recently, Gierz et al have also ob- 
served a large and anisotropic spin splitting in graphene 
samples on SiC substrate. 

The microscopic Hamiltonian describing the RSO cou- 
pling in graphene, shown in Eq. (jlbl) . was introduced by 
Kane and Mele in Ref. Il3i and its prope rties were subse- 
quently studied in e.g. Refs. [qlg jl^ila . Here we report 
on an interesting feature of the low energy spectrum of 
this Hamiltonian which has attracted little attention so 
far but which can be important in systems where RSO 
coupling is much stronger than ISO. Examples for such 
systems appear to be the ones studied in the experiments 
of Varykhalov et al^ and Gierz et al^. We show that close 
to the K and K' points of the graphene's Brillouin zone 



the Hamiltonian of monolayer graphene with RSO inter- 
action is in good approximation equivalent to the Hamil- 
tonian of bilayer graphene>12> and that RSO interaction in 
monolayer graphene leads to trigonal warping (TW) of 
the energy bands. This equivalence of the Hamiltonians 
may seem surprising since the two physical systems are 
very different. We show numerically that as long as RSO 
interaction is stronger than ISO, it affects the band struc- 
ture in two important ways: i) TW changes the topology 
of the low energy bands at the K (K') points which may 
be detected as a tree-fold increase in the minimal con- 
ductivity of the sample at very low electron densities, ii) 
RSO is also manifested through the non-isotropic spin- 
splitting of the bands and in particular, through the non- 
constant band splitting at the K (K') points. 

Tight binding model. — We start our discussion by con- 
sidering the tight binding (TB) Hamiltonian of mono- 
layer graphene which includes RSO coupling^ as well. 
It is given hy H — Hq + Hji , where 






(la) 



Here 70 is the hopping amplitude between nearest neigh- 
bor atoms, a]^ (ctio-) creates (annihilates) an electron in 
the ith unit cell with spin a on sublattice A, while b- 
(bja) has the same effect on sublattice B. The sum is 
taken over nearest neighbor atoms {i,j) and h.c. stands 
for hermitian conjugate. The Hamiltonian Hn describes 
the Rashba spin-orbit coupling and it readsi^: 
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where s = (sx^ Sy, Sz) are the Pauli matrices representing 
the electron spin operator and /i, i^ — 1,2 denote the fiv 
matrix elements of the Pauli matrices. Moreover, d^^ j^ — 
d(ij)/(i are unit vectors, where <i{i,j) points from atom j 
to its nearest neighbors i and d = |d^jj)|. The strength 
of the spin-orbit coupling is denoted by Xji which may 



arise due to a perpendicular electric field or interaction 
with the substrate. 

Using the operator aj^^ = ^« X^i '^"'^''^L' where R.; 
is the Bravais vector of the ith unit cell and q lies in the 
first Brillouin zone (and similarly, introducing (3^^ acting 
on sublattice B) it is easy to find that 
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Here ai , a2 are the lattice vectors of graphene and it was 
convenient to introduce the vector as = 0. Finally, the 
matrix representation of the Hamiltonian H = Hq + Hr 
acting on wave function ^ = (^Ati ^st' *A^, "^Bi) can 
be written as 
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where I?±(q) = ±Dx{q) - iDy{q). 
nian can be used for calculations in the whole Brillouin 
zone (BZ). As it has been shown in Ref. 14, the Hamil- 
tonian in Eq. ([3]) preserves the particle-hole symmetry 
therefore in the following we will focus on the conduction 
bands only. One finds that the electron and hole bands 
touch at the K and K' points of the Brillouin zone (here 
K = (2b2 + bi)/3 and K' = (2bi + b2)/3, where bi 
and b2 are the reciprocal lattice vectors of the graphene 
sheet). Next, it will prove instructive to derive a low en- 
ergy Hamiltonian which describes the excitations in the 
vicinity of the band-touching points. 

Continuous model. — We expand the Hamiltonian ^ 
around the Dirac point K of the BZ in terms of k == 
q — K for |k| w 0. Then one can perform a unitary 
transformation using the matrix U = (-^-^2 + ^^'^z) ^ 
Sz (here /2 is a 2 by 2 unit matrix) and take the inverse 
Fourier transformation corresponding to the replacement 
Kk —i' f) = —iKV (here p = {px,Py) is the momentum 
operator). One obtains finally the following Hamiltonian 
for monolayer graphene with RSO interaction: 



Hk^ 



vfP- vxp^ 

vfP+ SiXr 

3iAfl vpp_ 

\V\p- vfP+ 
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Here vp — 37o(i/(2?i), v\ — 3Afld/(2fi,) andp± = px^ipy 
Apart from a unitary transformation this Hamiltonian is 
the same as the Hamiltonian of bilayer graphene, if the 



asymmetry between the on-site energies of the bottom 
and top layers and the small hopping amplitude 74 can be 
neglectedi^. This conclusion may seem surprising since 
the TB models of the two systems look rather different. 
As vx <^ VF^ when calculating the band structure the ma- 
trix elements {Hk)^ and {Hk)^i are small compared to 
the matrix elements which are proportional to vp and one 
can question if it is important to keep them. Indeed, the 
theoretical treatment of RSO in Refs. I4ll8ll9lil2l is based 
on a Hamiltonian in which {Hk)^ = {Hk)^ = 0. Note 
however, that these terms cause trigonal warping of the 
bands at low energy and they can therefore be important 
since TW results in the change of the topology of the en- 
ergy bands near the Dirac point. This can be easily seen 
by recalling the band structure of bilayer graphene with 
TW^^ or by looking at the plane wave solutions of the 
Schrodinger equation ^^^^(k) = _B(k)*(k). The four 
eigenvalues of the Hamiltonian (|3]) as a function of the 
wave number k = A: (cos a, sin a)"'" are given by 



£;±(k) = ±hvp y 2 [fc2 + P (2 + /32) + (-1)«VT 

"(5a) 

where /3 — vx/vp — A/7 is the dimensionless strength of 
the spin-orbit coupling, kx — 2f3/d, n — 1,2, and 

T = fct + 2k^kl (2 - 13^) + k^p^iA + /3^) - Sk^kxfi sin(3a). 

(5b) 
Figures [Ua) and (b) show the contour plot of the con- 
duction bands E^Ck) and i?2^(k). At k = the split- 
ting of the bands is SXr [see Fig. Hfc)]. Similarly to 
bilayer graphene, the bands have a threefold rotational 
symmetry. For energies smaller than the Lifshitz energy 
El = 3Afl,32/(4-|-/32), the constant energy lines of i?+(k) 
are broken into four pockets, which can be referred to as 
one central and three leg parts [this is shown in the inset 
of Fig. [Ija)]. The eigenvalues ±E^{\i.) become zero at 
the K point of the Brillouin zone, ie, at fc = 0, and at 
the center of the three leg parts located at fc = kgo and 
a — — 7r/6;7r/2; 77r/6. The distance between the center 
of the leg parts and the Dirac-point (K point) in con- 
tinuous model is fcgo = /^^ f i which is equal to /ego in 
Ref. 14 for the experimentally relevant case of A/j <C 70- 
The authors of Ref. [ij have also noticed the three addi- 
tional zero-energy points around the K point but it was 
not recognized that the Hamiltonian ^ of monolayer 
graphene with RSO is equivalent to the Hamiltonian of 
bilayer graphene with TW. An important difference how- 
ever between bilayer graphene and monolayer graphene 
with RSO is that the splitting A£'(k) = ^2^(k) - S+(k) 
of the two conduction bands is much smaller in the latter. 
This energy is in the range of 0.3 — 0.48eV-'^ for bilayer 
whereas a recent experiment^ shows that the splitting is 
about 13meV for RSO in monolayer graphene. There- 
fore in the latter system even for small electron densi- 
ties states in both E^Ck) and i?^(k) would be occupied, 
whereas in the case of bilayer graphene the occupation of 
the upper band can usually be neglected. 

So far in our discussion we have completely neglected 




A5O=0 



AsO«^R 



2AsoS3Xr 2A,o>3V 



FIG. 1: (Color online) The k dependence and spin polariza- 
tion of the energy bands: (a) E^(k) and (b) _E^(k) around 
the K point for /3 = 0.034. Solid lines show the isoenergy 
contours, the arrows indicate the direction of the spin at cer- 
tain points in the BZ. The wave vector components kx, ky 
are in units of 3Afl/(2ft«F). The inset in Fig. [ija) shows a 
close-up of E^{k) in the vicinity of the K (K') point. The 
pocket structure has a 27r/3 rotational symmetry and the dis- 
tance of the pocket minima from the K (K') point is fcgo 
(see the text). The dashed line Fig. [ijb) is a circle, it in- 
dicates that the inosenergy lines of -Ej(k) are also slightly 
distorted, (c)-(f): a cut through the energy bands for kx — 
and (c) Aso = 0, (d) Aso < Xr, (e) 2Aso < 3Xr and (f) 
2Aso > 3\r. Solid blue lines denote electron bands, while 
dashed red lines correspond to valance bands. 



the intrinsic spin-orbit coupling ISO. As the Lifshitz en- 
ergy El itself is a very small energy scale, it is impor- 
tant to understand whether the trigonal warping and 
the pocket structure survive if a small, yet finite ISO 
is also present. To this end we have numerically cal- 
culated the band structure using the Hamiltonian H = 
Ho + Hr + Hso with Hso = ^soo-zSz (see Ref.Q) where 
az = 1(— 1) describes state on sublattice A (B) . The re- 
sults of the calculations are shown in Figs. [Ilc)-(f). For 
Xr ^ Ago the topology of the lower energy band around 
the band touching points is preserved [though the bands 
get a global shift, compare Fig. [TJc) and (d)]. Further- 
more, if 2A50 < 3A/i;, most of the characteristics of the 
RSO spectrum survive, but the minima of the leg pock- 
ets will be at slightly higher energy than the minimum of 
the central pocket [see Fig.[IJe)]. Finally, if the ISO split- 
ting become larger than SA/j, the pocket structure van- 
ishes completely [Fig.[T][f)]. Our findings in Figs.[TIc)-(d) 
are in good qualitative agreement with the density func- 
tional calculations of Ref. |j except that in the results of 
Ref. y the TW seems to be absent. This is not surpris- 
ing however, since for the weak Rashba splitting assumed 



there (3A_r w 25/xeV) trigonal warping is negligible. For 
stronger Xr values and if Xr » Ago, which is likely to be 
the case in the experiment of Refs. 00, TW is expected 
to be more important. 

We now turn to discuss some of the experimental impli- 
cations of our results and compare the predictions based 
on the Hamiltonian ^ and ^ to some of the results of 
recent experiments. One can immediately see that if the 
RSO coupling is much larger than the ISO, our results 
predict that due to the pocket structure shown in the in- 
set of Fig.[T]^a) for very low electron densities the minimal 
conductivity is three times larger than what one would 
obtain by neglecting TW. (For the relevant calculations 
for bilayer graphene see Refs. llTlIlSl ) 

Another way to study the TW experimentally would be 
to scan the spin splitting of the bands across the BZ. By 




FIG. 2: (a) The spin splitting along the TK line of the BZ 
using the Hamiltonian of Eq. ((3|. The inset shows the com- 
parison of the experimental data from Ref. and the calcu- 
lated band splitting with finite TW (solid line) and with no 
TW (constant dashed line), (b) The band splitting along 
TM for Xr, = 80meV, A^^ = 60meV, Xr^ = 40meV 
and Xr^ — lOmeV obtained from Hamiltonian ((Sj. Verti- 
cal dashed line indicates the minimum of the curves. The TB 
hopping parameter was 70 = 3eV. 

looking at Figs.[lja) and (b) it is easy to see that the split- 
ting is anisotropic around the K point and would reflect 
a threefold rotational symmetry. This is an irnportant 
difference compared to the modell used in Refs. HlsUolil^. 
In contrast to the measurements of RefT?, in this model 
the splitting does not go to zero along the TK line, which 
may indicate a strong influence of the SiC substrate that 
the theoretical models considered so far do not capture. 
Since the Hamiltonian ^ can only be applied around 
the K point, in order to explore a larger part of the BZ 
we also performed calculations using the TB model [see 



Eq. ([3])]. The band splitting along the TK line assum- 
ing RSO coupling constant SXr = 12meV is shown in 
Fig. [2ja). The most important feature to be observed is 
that the RSO splitting in the vicinity of the K point is not 
constant as predicted by the model which neglects TW, 
but has a small curvature. This can also be seen in the 
inset of Fig. ^a.) which shows the comparison of the cal- 
culations and of the data from Fig. 4(c) of Ref. i6|. The ex- 
perimental data are consistent with the calculated curve, 
unfortunately the measurement resolution does not allow 
to make a definitive judgment as to whether the splitting 
is constant or not. In Fig. ^h) we plot the spin splitting 
along the TM line in the BZ for various RSO coupling 
strengths. The maximum of the band splitting equals 
3Ai^ and it is always attained at the same k wavenum- 
ber, not depending on A/j. It is interesting to note that 
the shape of these curves is very reminiscent of the mea- 
surements shown in Fig. 3(c) of Ref. |2(/, which were made 
on graphenc on Ni substrate. The reason of the splitting 
in this system, however, appears to be uncleari2ii2^. 

The SARPRES techniquei^ can not only probe the 
band structure but also the spin polarization of the 
bands. The results of Ref. |a suggest that the spin polar- 
ization has circular symmetry around the K point and 
that the spin is polarized in clockwise (anticlockwise) di- 
rections in the two bands at the Fermi energy. Our cal- 
culations for the spin polarization of the energy bands 
EJ' (k) and E2 (k) are in agreement with this experimen- 
tal result [see in Figs.[T](a) and (b), respectively]. It might 
seem surprising that the spin polarization of the band 
E^ (k) shows a simple rotational symmetry, without any 
signs of TW observed in the spectrum. As we will demon- 
strate below, TW affects the spin polarization only in 
higher orders of k. The calculation of spin polarization 
goes along the lines of a similar calculation in Ref. jsl. The 
expectation value of the spin components on sublattice A 
is given by s° = h/2{sin9'^ cos ip'^, sin 9"^ siiup'^, cos 6 f)"^ . 
We find that in the lowest orders it depends on fc = |k| 



and 
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where the index n = 1,2 correspond to energy eigenvalues 
£:+(k) (or E-{k)) and £;^(k) (or £;2"(k)) respectively, 
and K — f3/d. Similar calculations for sublattice B yield 
6'^(k) = TT - 6'°,(k) and ip'^{k) = (/?°(k). As one can 
see from Eqs. (pa|) and ([Fb]) . TW does not play role in 
leading order of fc. This explains the circular symmetry 
of the spin polarization in Figs.[IJa) and (b). Moreover, 
the component of the spin perpendicular to the graphene 
sheet points into opposite directions on sublattices A and 
B. Therefore, as in Ref. @, we find that the average spin 
s = (s° -I- s^)/2 is polarized in the plane of the graphene 
sheet, it is perpendicular to the wave vector k and its 
magnitude grows up from zero (at k = 0) to h/2 further 
from the K point, see also in Fig. [Ha) and (b). 

In summary, we have studied the RSO coupling in 
graphene monolayer. We have found that in case both 
ISO and RSO interactions are present, as long as the 
RSO coupling is stronger, the low energy bands will show 
TW in the vicinity of the K point of the BZ. Our pre- 
dictions for the spin splitting of the bands along certain 
high symmetry directions in the BZ are consistent with 
the results of a recent experiment^ which measured the 
band splitting due to RSO interaction. Finally, we have 
demonstrated that the TW of the bands does not affect 
the spin polarization in the lowest order of the momen- 
tum therefore it shows rotational symmetry around the K 
point. These results might be relevant for other novel low 
dimensional systems where RSO coupling is importanti^. 
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